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We study how the coherence of noisy oscillations can be optimally enhanced by external locking. 
Basing on the condition of minimizing the phase diffusion constant, we find the optimal forcing 
explicitly in the limits of small and large noise, in dependence of phase sensitivity of the oscillator. 
We show that the form of the optimal force bifurcates with the noise intensity. In the limit of small 
noise, the results are compared with purely deterministic conditions of optimal locking. 
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Autonomous self-sustained oscillations may be ex¬ 
tremely regular (like, e.g., lasers) or rather incoherent 
(like many biological oscillators, e.g., ones responsible 
for cardiac or circadian rhythms). A usual way to im¬ 
prove the quality of oscillations is to lock (synchronize) 
them by an external pacing 00 - This is used in radio- 
controlled clocks and in cardiac pacemakers; also circa¬ 
dian rhythms are nearly perfectly locked by the 24-hours 
day/night force. 

In this letter we address a question: which periodic 
force ensures, via locking, the maximal coherence of a 
noisy self-sustained oscillator? Of course, one has to fix 
the amplitude of the force, so the nontrivial problem is 
in finding the optimal force profile. We will treat this 
problem in the phase approximation 0 ], which is valid 
for general oscillators, provided the noise and the forcing 
are small. In this approximation the dynamics of the 
phase reduces to a noisy Adler equation @,0 , and the 
maximal coherence is achieved if the diffusion constant of 
the phase is minimal. It should be noted that an optimal 
locking problem has been recently discussed for purely 
deterministic oscillations. There, the optimal condition 
was formulated as the maximal width of the Arnold’s 
tongue (the synchronization region) or as the maximal 
stability of the locked state @-0. In our case there is 
an additional parameter, the noise intensity, and we will 
show that the optimal force profile depends on the noise 
amplitude. Below we will also compare the limit of small 
noise with purely deterministic setups. 

Let us consider a self-sustained oscillator with fre¬ 
quency w, its phase in presence of a small Gaussian white 
noise obeys the Langevin equation 

^ w + p-mm) = m -(i) 

where /3” 1 is the noise intensity. A small periodic forcing 
with frequency f l leads, in the first order in the force 
amplitude, to the following phase dynamics 00 : 

w + s(tp)f(Qt) + P 1,/2 £(i)- (2) 

Here s(<p) is the phase sensitivity function (a.k.a. phase 
response curve), and fiflt ) is the phase-projected force 


term. Our goal will be to find such a forcing /(•) that 
maximizes the coherence, i.e. minimizes the diffusion 
constant of the phase ip. This optimal force will depend 
on the phase sensitivity function s(-) and on the noise 
intensity p. 

As the first step we introduce the slow phase (f> = ip — 
fit and perform the standard averaging over the period 
27 tH _ 1 00, this yields 

f = u - si+ m +r 1/2 «t) = -^ + r l/2 m 

( 3 ) 

where 

9 ^ = 'hJ 0 c M^+y)/(y)> ( 4 ) 

and we introduced the “potential” 

v(<t>) = {Q,-u)<t>- j g{y)dy. ( 5 ) 

Let us consider a situation, where the mean frequency 
of oscillations is exactly that of the forcing; this means 
that the slow phase <j> performs a random walk without 
a bias. This happens for a purely periodic, non-inclined 
potential. This condition, as it follows from ([5jl. defines 
the optimal frequency of the forcing 

Q = u + (s)(f), (6) 

where we denote (/(</>)} = (27r) -1 f(cj>) d(f>. Thus, 
without loss of generality we can assume that (s) = (/) = 
0 and fl = ui. 

The problem of finding the diffusion constant I? of a 
particle in a periodic potential v, driven by a white Gaus¬ 
sian noise, has been solved in Ref. [8] (and generalized to 
the case of an inclined potential in Ref. 0): 


(exp(/3u)}(exp(-^u))’ 

where Dq is the bare diffusion constant without potential. 
Thus, the problem of maximizing the coherence reduces 
to maximizing the expression 

C = (exp(/3f))(exp(—/?u)}. 


( 8 ) 
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As an additional condition we have to fix the intensity of 
the force: 


(/ 2 ) = const. (9) 

The formulated optimization problem is quite complex 
to be solved in general. Therefore, below we consider 
some simplifying cases, and will perform a rather full 
analysis for a simple bi-harmonic phase sensitivity func¬ 
tion. The main feature we will focus on, are bifurcations 
in dependence on the form of this function and on the 
noise intensity; we will see that different forcing wave¬ 
forms provide optimal coherence in different domains of 
the parameter space. 

For the analytical consideration below it is convenient 
to use Fourier transforms, which wc will denote by capi¬ 
tals: 

_ 1 t 27r 

s(x) = > Sk exp[ifc:r], Sk = — / s(x) exp[— ikx]dx, 

k 2nJ ° 

. ( 10) 

and the same for functions f,g,v, Fourier harmonics of 
which we denote as i4, Gfc, 14, respectively. Because g{cj)) 
is according to © a convolution of / and s, and v is the 
integral of g , we have 

Gk = SkF-k, 14 = ik~ 1 SkF-k■ ( 11 ) 

The condition on the norm of the force © now reads 

F |Ffc| 2 = const. (12) 

k 

We start with the case of strong noise (small /3)- Ex¬ 
panding © , we obtain a simple expression for the quan¬ 
tity to be maximized: 


Suppose now that v m i n = v(x 2 ) and v max = v(x\). Then 


In C ~ Vmax Vmin — 
r 2ir 


roc 2 

/ 9{<t>) 

J Xi 


d(f) = 


Xi 

2 r'Z 7T />X+X 2 

= 7T f(x)dx / s(y)dy. 
^ 7r Jo Jx+X-L 


(17) 


Using additionally condition © with a Lagrange multi¬ 
plier, we obtain 

t-X+X 2 

f{x) = const / s{y)dy. (18) 

J X-\-X\ 


Substituting this into conditions g{x i) = g(x 2 ) = 0, we 
get an equation for A = X 2 — x 1 (only this difference is 
important, but not the values of X\,X 2 )'- 


p ATT /»ZA 

p{ A) = / s(z) / s(z + y) dz dy = 
Jo Jo 

= 2nJ2k- 1 \S k \ 2 sinfcA = 0. 


(19) 


This equation has always a solution A = 7r, but depend¬ 
ing on the form of the phase sensitivity function s there 
can be other solutions, corresponding to local maxima of 
C; one has to compare different possible values of A to 
find the global maximum. Once A is found, the corre¬ 
sponding force can be expressed as 

r t ~& exp '^ A l- 1 . (20) 


Here below we present the simplest nontrivial example, 
where it is possible, in addition to the asymptotic cases 
of small and large noise considered above, to perform the 
analysis for intermediate noise levels. We consider the 
bi-harmonic phase sensitivity function 


C « 1 + P 2 (v 2 ) = l + p 2 J2 k~ 2 \S k \ 2 \F k \ 2 . (13) 

k 

Together with condition the maximum can be found 
by virtue of Lagrange multipliers: 

\Fk\~Sk,K, where K = arg max(fc _2 |5 fc | 2 ). (14) 

Thus, for large noise, the optimal forcing is purely har¬ 
monic one f(x) ~ cos(Kx), where K is determined 
from (fill) . 

The case of small noise is the limit /3 —> 00 . In this case 
the integrals in the expression © can be asymptotically 
estimated as Laplace integrals: 


s(x) = 2^/q cos x + 2y/l — geos 2x, (21) 

where parameter q describes the relative weight of the 
harmonics. 

The limit of strong noise flTjJ) , with |Si| 2 = q, IS 2 I 2 = 
1 — q, yields 


f(x) 


cos 2x if 0 < q < 1/5, 
cosa; if 1/5 < q < 1. 


( 22 ) 


The limit of weak noise leads to the following expres¬ 
sion for function ©: 

p(A) = g sin A -|-^—^sin2A. (23) 


(exp(@v)) « (2tt) 1 exp(/3vmax), 
(exp(-/3u)) ~ (2 tt) -1 exp(-Pvmin), 


what gives 


C exp(f3(Vmax XmlnY}- 


(15) 


For q > 1/2, the only root in C3l) is A = n, while for q < 
1/2 there is an additional root Ai = arccos(— q/(l — q)). 
Substituting this into (0(71) . we obtain for small noise 


(16) 


= 1 -= I 2 " ,f0 f' I<1/2 ’ 

11 1 21 \l if 1/2 <g< 1. 


( 24 ) 
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FIG. 2. (color online) The intensity of the first harmonic 
component in the optimal force as function of parameter q, 
for three optimization criteria: top solid red line: maximal 
coherence in the weak noise limit El; middle dashed green 
line: maximal width of the synchronization region OU1) : bot¬ 
tom dotted blue curve: maximal linear stability of the locked 
state 


Unfortunately, after substitution of this potential in the 
expression J5J for the factor C, we obtain integrals which 
cannot be expressed in a closed analytic form. However, 
for a purely first-harmonic forcing (6 = c = 0 ) and a 
purely second-harmonic forcing (a = 0), the factor C as 
well as its derivatives can be expressed via first order 
Bessel functions. Thus, it is possible to find the domains 
of stability of these pure forcing terms analytically, for 
arbitrary values of noise intensity /3. These lengthy but 
straightforward calculations give the stability boundaries 
in a parametric form: The first-harmonic force loses sta¬ 
bility at the curve on the (/3, q) plane given according 
to 


q = 


z(-h{z) + Io{z)) r = — 

8Ii(z) + z(-h(z) +Iq(z)) ’ y/q' 


(27) 


FIG. 1. (color online) (a) Domains on the plane of parame¬ 
ters (q, P) where the optimal force has one harmonics (2: the 
second one, 1: the first one), and two harmonics (1+2), ac¬ 
cording top expressions H27I28I) . (b) The intensity of the first 
harmonics |Ui | 2 as a function of q for different noise intensities 
P. Thin dashed black line shows the limit of small noise El- 
(c) Two dependencies of |Fi| 2 on the noise intensity P showing 
bifurcations from one-mode to two-mode solutions at critical 
values of p. 


Let us now consider general noise intensities. The 
forcing in this case should be also generally bi-harmonic 
(higher harmonics disappear according to (11 ID 1: 


f(x) = acosx + b cos 2x + c sin 2x, 


(25) 


with unknown constants a, 6 , c satisfying a 2 + b 2 + c 2 = 1 . 
In this representation the potential v(x) reads 

/ \ r- . V 1 ~ Q b ■ o V 1 - qb o 

v(x) = —y/qasmx ---sm 2 x---cos Ax. 

( 26 ) 


The stability boundary of the second-harmonic solution 
is 

= h(z) = 2z 

h{z) +2zl 0 (z) ’ y/1 — q 

We illustrate these domains in Fig. [L] Here we also show 
numerically obtained dependencies of |Fi | 2 (the intensity 
of the second harmonis is I-F 2 I 2 = 1 —|Fi| 2 ) on parameters 
q and /?, demonstrating bifurcations on the form of the 
forcing. 

Next we discuss a relation between different criteria 
used for the “optimal locking”. While here we optimize 
the coherence in the presence of noise, in Refs. [1, H 
purely deterministic criteria have been suggested. It is 
instructive to compare them with our approach in the 
limit of small noise. Suppose that the coupling function 
g{(j)) has zeros at d +2 (where </>i is the stable one), and 
extrema at ^> 3 , 4 . In the approach of Q, the linear stabil¬ 
ity at the stable equilibrium |g , (aq)| is maximized. In the 
approach of |J, the width of the synchronization region 
~ \g(4>3) — g(4> 4 )| is maximized. In our maximization of 
the coherence, the potential barrier for a noise-induced 
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phase slip 


Jj* £Kz)cte 


should be maximal. For the 


discussed above example of a bi-harmonic phase sensitiv¬ 
ity function HD, all the optimal forcings can be found 
analytically, they are generally also bi-harmonic. The 
approach of fij] yields in this case 


|Fi| 1 2 3 4 5 = l-|F 2 | 2 = g/(4-3g), 


while the approach of 3 gives 


|Fi| 2 = 1 - |F 2 | 2 


3^ if 0 < q < 4/5, 
1 if 4/5 < g < 1. 


(29) 


(30) 


In conclusion, we have studied the problem of maximiz¬ 
ing coherence of oscillations by external locking, in the 
phase approximation. The optimal phase forcing func¬ 
tion depends not only on the phase sensitivity function of 
the system, but also on the noise intensity. For large noise 
a purely harmonic forcing is optimal, the number of the 
harmonic depends on the phase sensitivity. For smaller 
noise, a bifurcation to a more general, multi-harmonic 
forcing may occur. We have also demonstrated, that dif¬ 
ferent optimality conditions in the purely deterministic 
case lead to different optimal forcing functions, which 
also differ from the limit of small noise when optimiza¬ 
tion of the coherence is performed. 


We compare the results in Fig. [2] One can see that for 
the minimal coherence, presence of a strong first harmon¬ 
ics component in the forcing is more important than for 
other criteria. 
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